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A surface ship subjected to an underwater explosion is exposed to shock waves 
over a short period of time which can vary in magnitude based on charge type, size, and 
location. The energy of those waves impinging upon the hull is transmitted throughout 
the ship’s structure and vital equipment. The dynamics of the shock waves also influence 
the fluid surrounding the outer hull of the ship, creating an area of cavitating fluid. The 
combination of the shock waves, bubble pulsations, and cavitating fluid induce shipwide 
vibrations on hull supports and mission essential equipment which may become 
inoperative. In view of congressional requirements for new ship designs and systems to 
be shock tested, this thesis investigates the modeling of a preliminary design (Flight I) of 
the Arleigh Burke Destroyer (DDG 51) exposed to an underwater explosion. The effects 
of cavitation on one and two dimensional models is explored to determine if cavitation 
effects are substantially important to a three dimensional ship model. Validation of 
modeling underwater explosion effects upon a ship model can provide potential insight 
and savings in cost for future live fire testing and evaluation of the Flight IIA (DDG 79) 
design of the Arleigh Burke Destroyer. 
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I. 



INTRODUCTION 



A. BACKGROUND 

Today’s U.S Navy consisting of technologically advanced ships such as the 
Ticonderoga class cruiser and the Arleigh Burke class destroyer are being driven to 
operate in the littoral, close to shore, environment as dilineated in the Department of the 
Navy’s strategic planning document [Ref. 1]. These platforms equipped with the Aegis 
weapon system were designed for open ocean, “blue water” tactical scenarios but due to 
the end of the cold war, a shift in emphasis is placed on low intensity conflicts which 
occur near the sea-land interface. This has increased the potential for “cheap kill” 
opportunities against U.S. warships such as an underwater mine as ships close the shore 
and the availability of such weapons is high and the cost low for third world country to 
purchase. 

As recently as 1991 during Operation Desert Storm, the USS Princeton (CG-59), a 
Ticonderoga class guided missile cruiser, struck a floating mine near the bow. Although 
the rupture to the hull was localized, the subsequent whipping of the hull caused severe 
hull girder damage near the stem of the ship and weapon systems to go off line. This 
type of tactical environment now places a greater responsibility on ship designers and 
engineers to meet requirements as specified in MIL-S-901D [Ref. 2], 

A number of studies have been conducted on submerged cylinder type models 
exposed to an underwater explosion in order to predict the hull response of structures 



such as a submerged submarine. These studies did not include the effects of cavitation on 
the wet surface of the model or of the cavitation which occurs at the free surface. For 
surface ship type studies, simplified beam models such as that conducted by Hicks [Ref. 
3] have concentrated only on the effects of whipping of surface vessels due to gas 
bubble dynamics. In this age of improving computer hardware, larger memory and 
software programs such as the USA (Underwater Shock Analysis) code with the CFA 
(Cavitation Fluid Analyzer) [Ref. 4] to model underwater shock effects and improving 
CAD programs with built in FEM (Finite Element Model) generators such as IDEAS 
[Ref. 5] and MSC/NASTRAN [Ref. 6], the engineer has a greater ability to improve ship 
modeling and predict the vibrational response of ship hulls and internal equipment. 

B. SCOPE OF RESEARCH 

This thesis studies the methods of modeling one and two dimensional finite element 
models with and without an external fluid. Two different computational methods are 
selected in the USA code for each model arrangement in addition to incorporating the 
CFA module in the solution process. The goal is to determine the cavitation effects on a 
surface vessel and develop the most simplest and conservative method in modeling the 
underwater explosion. A one dimensional model is created from Bleich’s paper [Ref. 7] 
which defines the conditions of a square plate sitting on a fluid surface and a fluid 
column. In additon, an exponetial decaying plane wave models the underwater explosion. 
A second finite element model is created to represent the cross-sectional area of a simple 
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ship hull. This two dimensional model includes the surrounding fluid to the waterline of 
the hull and introduces the effects of a free surface due to an actual charge placed at a 
distance from the hull. The results from these first two models provides some 
simplification for the three dimensional notional hull structure of the Flight I of the 
Arleigh Burke. Although live fire testing and evaluation has already been conducted on a 
Flight I design, a Flight IIA design starting with DDG-79 possesses hull differences 
shown in Fig.l which may require shock testing. Improvements in the area of modeling a 
surface ship such as the Flight IIA design from existing CAD programs and fluid 
behavior during an underwater explosion through the USA code, can create potential 
savings in future Flight IIA shock trial tests. 




Figure 1. Major Changes (Flight II to Flight IIA) [Ref. 8] 
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II. UNDERWATER EXPLOSION 



A. FLUID BEHAVIOR 

When an underwater explosion occurs, the solid explosive material (such as 
HBX-1, TNT, etc.) suddenly reacts, leaving behind high temperature and pressure of 
gaseous products. An initial wave of compression is produced due to a sudden 
discontinuity of pressure which can reach a peak pressure of 2 x 10 6 psi. This steep 
fronted wave known as a “shock wave” propagates radially outward at a speed several 
times the speed of sound in water and roughly decays exponentially. The velocity of the 
shock wave decreases such that at approximately 10 times the explosive charge radius, 
the disturbance is essentially traveling at the acoustic velocity in water, 5000 ft/sec. [Ref 
9-.pp 3-7 ] 

The pressure sensed at a specific point in the fluid or at a structure is essentially a 
step increase to the peak pressure followed by an exponential decay with time until the 
pressure has dropped to approximately one third of the peak pressure. After falling to 
approximately one-third of the peak pressure, the pressure level falls off inversely with 
distance [Ref 9:pp3-7], 

The cavity of gaseous products left behind at high pressure forms a bubble which 
subsequently expands to relieve the difference in pressure, accelerating the surrounding 
fluid particles. The bubble continues to expand beyond the point of hydrostatic 
equilibrium (due to the inertia of the surrounding fluid) until a point of dynamic 
equilibrium is reached. The bubble then reverses its motion, continuing to contract until 
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dynamic equilibrium is again reached, where it quickly rebounds and again begins to 
expand as illustrated in Fig. 2. This oscillating bubble expansion and contraction 
continues until the energy of the reaction is fully dissipated or the bubble finally reaches 
the surface, venting the gaseos by products of the explosion. As the bubble rebounds, it 
greatly accelerates the surrounding water, generating a substantial pressure pulse (known 
as the bubble pulse). This bubble pulse can impart significant loads on structures in the 
vicinity. 

In addition to the initial shock wave and bubble pulse, based on the location of the 
expolsive charge with respect to the surface and sea floor, a surface vessel can also 
experience a combination of the following pressure waves as illustrated in Fig. 3: 

1 . Free-surface reflection 

2. Bottom reflection 

3. Bottom refraction (not shown) 

The directions of free-surface and bottom reflection waves can be determined through an 
application of Snell’s law at the surface and sea bottom, respectively. At the surface of a 
fluid, a reflected wave of negative pressure known as the “rarefraction wave” is formed 
with a value such that the sum of the direct and reflected pressures is zero along the 
boundary condition between air and water defined by the surface. This rarefraction wave 
travels through the fluid region shortly after the incident shock wave at any point which 
causes the incident shock wave pressure profile to be truncated at a point in time called 
“surface cut-off.” Bottom reflected and refracted waves are dependent on the 
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characteristics of the sea bottom and can contribute significant pressure waves in shallow 
water environments but are not considered for this study. 




Figure 2. Gas Bubble behavior 




Figure 3. Pressure Wave Profiles 
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Figure 4. Combined Pressure Profile of an Underwater Explosion 

The oscillating bubble pulse can also introduce a low-frequency forcing function that can 
resonate a ship girder frequency especially when the period of the bubble pulse is in the 
vicinity of the natural frequencies of the ship girder. Fig. 4 illustrates a combination of 
the pressure waves. 

A number of empirical relationships have been derived through numerous 
underwater explosion tests since World War II which provide a means to calculate 
explosion characteristics such as incident wave pressure (P max ) which impinges upon a 
target, shock wave decay constant (0), pressure as a function of time (P(t)), and bubble 
period (T) with derived constants for each respective explosive material provided: 
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t-t, 



P(t) = P„a«e 6 



1/3 



w 

P m ax=K,(— ) 



w 1/3 

0 = K 2 W 1/3 (—) A2 

R 



(psi) 

(psi) 

(msec) 



W ,/3 

5 (D + 33) 5 ' 6 (sec) 

W = Charge weight (lb) 

R = Standoff distance (ft) 

D = Charge Depth (ft) 
tj = Arrival time of shock wave (sec) 
t = Time of interest (sec) 
K 1 ,K 2 ,K 5 ,A 1 ,A 2 = Shock wave parameters 



B. STRUCTURAL BEHAVIOR 

The dynamic response of a linear-elastic structure can be defined, based on the 
number of elements and nodes that define the structure, as 

[M]{x(t)}+[C]{x(t)}+[K]{x(t)}= {F(t)} (2.1) 
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where [M] = structural mass matrix, [C]=damping matrix, [K]=stiffness matrix, 
{F(t)}=applied external force, {x(t)}= general displacement vector and derivatives with 
respect to time. For a structure in a fluid medium, the applied external force created by an 
underwater explosion, is defined as 

-[G ] [A fluid ] ({Pi(t)} + {P s (t)})= {F(t)> (2.2) 

where [G] = the transformation matrix between structural and fluid nodes, [A fluid ]= 
diagonal area matrix for the elements in the fluid, {P,(t)}= incident wave nodal pressure 
vector and {P s (t)} = scattered wave nodal pressure vector. 

The USA code provides the means to model the pressures involved in an 
underwater shock as a “forcing function” on a fully or partially submerged model. The 
code also utilizes variations of the DAA (Doubly Asymptotic Approximation) method 
[Ref. 10] that models the surrounding acoustic medium as a membrane on the surface of 
the structure actually in contact with the homogeneous fluid. The effects of cavitation on 
a structure modeled with a surrounding fluid can be also included through the CFA. The 
fluid motion and structural response is linked together through a series of compatibility 
relations that are matrix differential equations that relate the symmetric fluid mass matrix 
of the model’s wet surface, fluid particle velocities from the scattered wave normal to the 
structures wet surface, and rates of changing scattered fluid pressures. The USA code 
creates this matrix by utilizing an infinite, inviscid, incompressible fluid in an irrotational 
flow created by the structure’s wet surface and then creates a virtual mass matrix which is 
comprised of the added mass of the fluid on the elements of the structures mass matrix 
which are wet. Kinematic compatibility is obtained between the fluid motion and 
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structural motion by an invariance of virtual work and transformation matrices that 
create the final interaction equations for fluid and structure motion [Ref. 10]. 

Hence, the final equations of structure and changing fluid pressures create a set of 
“ augmented interaction equations” which the USA code solves, given a specified time 
increment, in an unconditionally stable manner. 



C. CAVITATION 



Two types of cavitation which can occur during an underwater explosion is that of 
“local cavitation” which occurs at the fluid-surface interface and “bulk cavitation” which 
occurs near the free surface that can extend beyond the target of interest by hundreds of 
feet. Utilizing Taylor flat plate theory, the particular models in this thesis will utilize only 
air backed plate conditions as illustrated in Fig. 5. 



I 



I 



Backed^; 



A i r 






I 



Pi(t) 



Vi(t) 






A i r 
Backed! 



Pr(t) 



Vr(t) 



Shock Wave Approaching 



Shock Wave Reflecting 



Figure 5. Taylor Plate Subjected to Plane Wave. 
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1 . 



Local Cavitation 



When the pressure pulse from an underwater explosion impinges upon a flexible 
surface such as the hull of surface ship, a fluid-structure interaction occurs which 
dynamically excites the ship structure. Looking at an infinite, air backed plate of mass, 
per unit area (m) subjected to an incident plane shock wave Pj(t), a reflection wave of 
pressure P r (t) will come off the plate. Letting v p (t) be the velocity of the plate and 
applying Newton’s second law of motion: 

mX TT= P r + P i (2.3) 

Letting the fluid particle velocities behind the incident and reflected shock waves be V;(t) 
and v r (t), respectively, the velocity of the plate becomes 

Vp(t) = V J (t) — v r (t) (2.4) 

Incident and reflected shock wave pressures are defined, respectively : 

PrpCvj (2.5) 

P = pCv r (2.6) 

where p=fluid density and C=acoustic velocity. Substituting the pressures (2.5) and(2.6) 
into (2.4) and utilizing for a high explosive shock wave 

__t_ 

p i(t) = p m,*e 9 (2.7) 
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where t = time starting after the shock wave arrives at the target, one can solve for P r : 



P, = ^ -pCv p = P max e e - pCv p 

Returning to the basic equation of motion (eq. 2.1) and substituting for Pj gives the 
following first order linear differential equation: 



dv 

m h pCv = 2P m e 

dt p 



(2.9) 



The above differential equation is solved for v p , 



V P = 



2P max e 



m(l-p) 



<-f) ( T> 

e u — e 0 



(2.10) 



where p=pC0/m and t>0. Total pressure on the plate is also found to be 



P + P = P 

A i x r x max 



2 -~i .-f 



e ° - 



1-P 1-P 



( 2 . 11 ) 



Hence, as P becomes larger, representing a light weight plate becoming more 
flexible, the total pressure will become negative at a very early time and as the total 
pressure approaches vapor pressure, local cavitation occurs in front of the plate. At this 
point, the plate is separating from the fluid and the maximum velocity for the plate is 
attained [Ref. 11] 
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2 . 



Bulk Cavitation 



The combination of the rareffaction wave and incident shock wave near the free 
surface creates a region of decreasing fluid pressure that approaches the vapor pressure of 
water. The rareffaction wave direction is determined by an image charge located an 
equivalent distance from the surface boundary in the air as the actual charge in the fluid. 
An upper and lower boundary begins to form which defines the bulk cavitation zone (Fig. 
6). The velocity of the water particle near the free surface will have primarily a vertical 
velocity at cavitation. As the reflected shock wave passes, the fluid is then acted upon by 
gravity and atmospheric pressure. 




Figure 6. Bulk Cavitation Zone [Ref. 11] 
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The upper cavitation boundary of Fig. 6 is the set of location points where the rarefraction 



wave arrives and the absolute pressure drops to the cavitation pressure. The total pressure 



which determines the upper boundary is the summation atmospheric, hydrostatic, incident 



shock, and rarefraction which can be defined as 



F(x,y) = K,( ) A, e~ ce +P A +Yy-K,( 

r, 



W U \ A 




r i=^(D-y ) 2 +* 2 



h =yj(D + y ) 2 +x 2 



where 



x,y = horizontal range and vertical depth of point 



r t = standoff distance from charge to point 
r 2 = standoff distance from image charge to point 
C = speed of sound in water 
D = charge depth 
0 = decay constant 

P A = atmospheric pressure 
y = weight density of water 
W = charge weight 



KjjAj = shock wave parameters 



The lower cavitation boundary is determined by equating the decay rate of the 



breaking pressure which is the rarefaction wave pressure that reduces a particular location 



of fluid to the point of cavitation pressure, approximately zero, to the decay rate of the 
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total absolute pressure of the fluid as illustrated in Fig. 6. This boundary can be defined 
as [Ref. 12] 



G(x, y) - - ~ 



Pl 

ce 



r 2 -2D(^) 



1 + 



^£l-a 2 -i 

. r i 





A, Pj 




f D + yV 




> 2 


r 2 - 2D 




J 


h 




l r 2 )_ 



+ 



^ D + y^ 
r, ) 



+-J-(P i+ P.+ry)=o 



where 



P i = K 



R J 



0 = K 2 W 1/3 



^W 1/3 ^| 

CkTj 



a 2 



Hence, any point that satisfies F(x,y)=0 and G(x,y)=0 determines the upper and 
lower boundary points. Fig. 7 and Appendix A provides a MATLAB code [Ref. 13] for 
bulk cavitation zones for various depths and charge weights for TNT explosive material. 
This particular tool was useful in modeling the 2D cross sectional hull and placing the 
charge in location that would minimize the effects of bulk cavitation and aid in 
concentrating effects of local cavitation. Fig. 7 illustrates the region of the bulk 
cavitation zone above the charge does not vary much with changing depth. However, the 
bulk cavitation zone increases in depth with appreciable increases in charge size. 30 feet 
below the keel (total depth of 50 feet) was used for the two dimensional study. 
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Depth(FT) Depth(FT) Depth(FT) 



100LBS TNT 






Figure 7. Bulk Cavitation Zones for a TNT charge located at the following depths 

—50ft, -100ft, -.150ft 
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III. COMPUTATIONAL FLOW 



Model preparation, translation, calculation, and post-processing involved a 
number of programs which is depicted in Fig. 8. Each process involves a number of 
output files which are input to the next step in the process or provides the user a means to 
verify results from each step before proceeding. 



MODELING 



TRANSLATION 



USA SOLUTION 



TRANSLATION 



POST-PROCESS 




Figure 8. Computer Program Flow 
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A. 



MODELING 



1. ASSET 

ASSET/MONOSC (Advanced Surface Ship Evaluation Tool/Monohull Surface 
Combatant) [Ref. 14] program is an interactive program which created the original ship 
dimensions for the three dimensional study. It provides a naval architectural tool to 
address ship design, which includes geometric definition of the hull and superstructure, 
hull subdivision, hull structure, resistance, appendages, equipment, and costs [Ref. 14]. 
An early design of the DDG 51 class destroyer with general weights of equipment 
groups and hull structure requirements based on inputs from designer imposed needs or 
mission requirements was exported out of ASSET in an IGES ( Initial Graphics Exchange 
System) file format which included general ship offsets which can be imported into 
IDEAS through the IGES translator. Although general bulkhead, deck, superstructure 
design was provided, no further details to subcompartmentized spaces, specific weights 
and locations of internal equipment, sonar dome offsets, and any other explicit details of 
the notional ship design was available. Appendix B provides a summary of the DDG 51 
notional model. 

2. IDEAS 

IDEAS, Master Series 2.0, [Ref. 5] was utilized to model all three cases of study. 
The IDEAS program is comprised of a large number of applications such as design, 
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modeling, testing, manufacturing, and simulation, each with a large group of subtasks for 
the user to utilize in icon panel form or menu driven selection. The program is also a 
“team concept” software which allows a number of users in a design or analysis team to 
take individual parts of the whole model to be designed, modified, or analyzed. For each 
model, the Simulation Application and all its respective subtasks were utilized to create 
(e.g. surfacing) and mesh the finite element model for export. 

a. One and Two Dimensional Models 

For the one dimensional and two dimensional cavitation study, one model 
was required to be created which included the structure and surrounding fluid. Fluid mesh 
thickness could not exceed the following requirement due to numerical stability of the 
CFA code: 

2 P~ D <5 

P.J., 

where p w ^density of water, p s =density of steel, t s =thickness of steel plate, and D= width 
of fluid elements. For example, utilizing p s /p w « 8, if the thickness of the plate or hull t s 
=0.5 in, D < 10 in. After creating all structure and fluid model regions, the following 
order for meshing was conducted due to the order that the CFATRAN and CFAGAL 
program organizes particular fluid nodes and elements: 

1 . Fluid volume elements (brick elements) 

2. DAA boundary (thin shell elements) 

3. Hydrostatic boundary (thin shell elements) 
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4. Specified Pressure regions (not utilized for this study, thin shell elements) 

5. Wet structural surfaces (thin shell elements) 

6. Remaining structural surfaces 

The thin shell elements for the DAA boundary, hydrostatic boundary, as well as for the 
specified pressure regions if used, provide a means to select the nodes in these particular 
regions that the CFA requires for numerical solution. The thin shell elements do not 
become a part of the soluiton process but act as “markers” in the CFATRAN program for 
the boundary element nodes which are specifically associated by the MAT field in the 
NASTRAN ASCII input file. Because of the numerical processing in the USA code, 
node numbering needs to be described in a counter-clockwise manner looking from the 
fluid towards the structure. Color coding particular surfaces and their respective elements 
or by creating specific groups of surfaces or elements in the meshing subtask-grouping 
icon greatly aids in changing particular characteristics such as the node numbering. 

A one psi pressure vector outward from the structure into the fluid is 
applied to all structural wet elements to act as a marker for the USA code fluid matrix 
development and restraint set conditions set for symmetric boundaries are created in the 
BOUNDARY CONDITION subtask. 

After meshing, two copies of the global mesh are made. One FEM model 
will be structure only with the one psi pressure, deleting fluid elements and nodes and 
one FEM model will be the fluid volume elements with the DAA, hydrostatic, and wet 
structure thin shell elements, deleting the other dry structural nodes and elements. These 
two FEM models are then exported out of IDEAS and import to MSC/NASTRAN. 
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b. 



Three Dimensional Ship Model 



Prior to the IGES translator task, modifications needed to be made to the 
actual IGES ASCII file for the correct unit system, ensuring IDEAS was preset to the 
same units (feet, inches, meters, etc.). 

Offset points and lines .were utilized as “wire guides” to recreate the model 
in IDEAS for creating all surfaces through lofting, mesh of curves, or planar creation 
shown in Fig. 9. The following sequence proved to be the best option in creating a 
complete three dimensional ship model in IDEAS, Simulation, Master surfacing task: 

1. HULL 

2. BULKHEADS 

3. DECKS 

4. SUPERSTRUCTURE: 

After all the surfaces have been created, all the edges must be connected or 
STITCHED in order for meshing to be successful and a check for free edges is performed 
in preparation for meshing. Any new surfaces that need to be created, such as internal 
bulkheads or connecting the superstructure to the hull, can be done through a surface 
bound by a combination of edges and three dimensional lines connecting those opposing 
edges. 

The required thin shell element properties are applied and a one psi 
pressure vector must be applied to all the wet surfaces through the BOUNDARY 
CONDITION subtask as well as any symmetry condition through the RESTRAINT SET 
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icon. IDEAS provides a full series of element checks (warping, distortion, interior angles, 
and normal directions) which must be done for proper solution processing. Only one 
FEM model was required for export since the CFA program was not applied to this 
particular model. 
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Figure 9. Imported ASSET DDG 51 Model 
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B. TRANSLATION 



Each of the following steps allows for output files to be viewed, verifying that 
information (number of nodes, eigenvalues, etc.) is being properly computed and created 
in the unformatted fortran files for the next step. 

1. MSC/NASTRAN 67.5 

The MSC/NASTRAN structure database is the required format for this particular 
version of the USA code during the initial start and restart. The USA code itself is 
conducting the solution process which would normally be a SOL 109 transient analysis 
process in MSC/NASTRAN. Appendix B provides the DMAP alters for IDEAS post 
processing which is attached to the beginning of the exported file out of IDEAS and for 
the MSC/NASTRAN restart, deleting unnecessary and unneeded lines up the bulk data 
section. This process generates the geometry information and the structure’s M, C, and K 
matrices 



2. CFATRAN/CFAGAL 

CFATRAN is a program modified from [Ref. 15] that takes the 
MSC/NASTRAN file of the fluid mesh model and groups the elements, nodes, and nodal 
connectivity in the required format for the CFAGAL program. CFAGAL will read the 
CFATRAN output and create a global access library of the fluid model and its 



26 



boundaries, just as the program FORGAL does for the structure model. This library will 
then be accessed during the FLUMAS process. 

3. NASSKY/FORGAL 

Nassky prepares the NAS TRAN UT1 output file which contains the geometry and 
matrix data of the structure into “skyline” matrix format , a generic Fortran unformatted 
file. FORGAL then converts this generic file data into a GAL file (Global Access 
Library) database which the USA code will access for structure related data. 

C. UNDERWATER SHOCK APPLICATION AND SOLUTION 

Appendix B provides sample inputs for the following programs which were 
utilized for the two-dimensional study. 

1. FLUMAS 

This is the first step in the USA code which creates the fluid mass mastrix of a 
structure wet surface in an infinite, inviscid, incompressible fluid by the boundary 
element technique. The key inputs here are the following [Ref. 16,17]: 

1. Mesh Geometry: fluid wet-surface and structure wet-surface mesh 

2. Element definitions: general curved surface or surface of revolution 

3. Material Property: mass density, speed of sound, DAA2 parameter 

4. Constraint options: location of free surface, half model, quarter model, long 
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cylinder, node reassignment in fluid-structure transformation if needed. 

For cavitation studies, the CFA is activated and the fluid volume database created 
by CFAGAL is incorporated with the FLUMAS output. Care should be exercised here 
during CFA solutions. The NCAV field refers to the number of elements on the wet 
structure plus the DAA boundary. During non-CFA runs, the NCAV field refers only to 
the wet elements on the structure since this defines the DAA boundary. 

2. AUGMAT 

This process in the USA code utilizes the equations specified earlier as the 
“augmented equations” by utilizing the data in the GAL file and FLUMAS output to 
construct the specific matrices. The key inputs are the following [Ref. 16,18]. 

1 . Fluid and Structure mass matrices 

2. Structural DOF correspondance table 

3. Fluid mesh geometry 

4. Fluid/Structure DOF transformation coefficients 

5. Type of solution to be conducted: plane wave, DAA2 options 

3. TIMINT 

This step constitutes the most time consuming step depending upon on the 
number of nodes, elements, and requested length of time to be studied. The 
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unconditionally stable staggered soluiton technique is utilized which basically involves 
the following [Ref. 16,19]: 

1 . Estimating the unknown structural restoring force [KJ (x(t)} at t + A t from the 
extrapolation of current and passed values. 

2. Solve the fluid equation and obtain a preliminary estimate of the total pressure 
vector at t + A t. 

3. Solve the structural equation for the displacement and velocity vectors at t + A t. 
Other steps are involved between these particular steps for matrix transformations and 
improved accuracy. General inputs from the user involve the following: 

1. Incident wave characteristics: location and type of charge source, standoff, 
specified pressure profile, exponentially decaying wave. 

2. Time step information: start and finish times, increments. 

3. Display data of particular nodes for checking output results in the solution 
file:desired node number and freedoms for displacement, velocities, or 
pressures on the DAA boundary. 

Although unconditionally stable, the user should utilize small enough time steps to 
capture accurate model behavior. Solution time history responses of all structural nodes 
and fluid nodes for all runs are stored in the file designated POSNAM in the TIMINT 
input. 
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D. POSTPROCESSING 



Two methods of post processing was done due to the growing amounts of 
required memory for larger models such as the three dimensional ship hull. 

1. XUDVT 

This second interface program converts the USA solution history in the 
corresponding MSC/NASTRAN data block, UDVT, a file format readable by 
MSC/NASTRAN. The initial MSC/NASTRAN job is restarted at the point just before the 
NASSKY translation, allowing MSC/NASTRAN to continue to a normal conclusion. The 
USA solution replaces that which would normally be generated by MSC/NASTRAN in a 
transient analysis. 

During CFA studies, the number of wet points required to be entered are those 
wet points totals generated by NASSKY which can be found in the nassky.out file. 
However, during non-CFA studies, the wet points refer to the total number of wet 
element faces. If the incorrect number of wet points or wet elements are entered, 
translation will fail or data history for some nodes will not be translated. XUDVT also 
allows displacement, velocity, acceleration, and pressure histories to be extracted in 
ASCII table format and placed in the PLT file for X-Y plotting by other programs (e.g. 
MATLAB). 
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2 . 



POSTPR 



This is a USA postprocessing program which allows extraction of particular node 
displacement, velocity, or pressure histories. Although not needed for non-CFA studies, 
this is the only current method which one can extract pressure history profile of fluid 
nodes during a CFA study. The data is placed in ASCII column tables which can be 
plotted by other programs such as Matlab. Due to the large solution files generated by the 
three dimensional ship model and the growing size of files during the XUDVT- 
MSC/NASTRAN-IDEAS translation and final post-processing, POSTPR proved to be a 
memory saving tool. The user must ensure the original UT1 file name is either deleted or 
changed in order to prevent a XUDVT failure since this process produces a UT1 file for 
the MSC/NASTRAN restart. 

3. MSC/NASTRAN (restart) 

Appendix B illustrates the required fields needed to restart MSC/NASTRAN with 
a new UT1 file created by the XUDVT process. This process utlizes the original 
model.MASTER file originally created by the initial MSC/NASTRAN run and creates a 
solution file for IDEAS to import (model.OP2) which contains all or selected time 
histories of displacement, velocity, acceleration, and stress. The user should utilize a 
different .OP2 file name from that of the .OP2 file created in the first MSC/NASTRAN 
run. Files usually triple in size compared to the POSNAM file output from TIMINT. 
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4. 



IDEAS 



The final .0P2 file can be imported into IDEAS through the FILE IMPORT- 
MSC/NASTRAN process. IDEAS will create a universal file (.unv) prior to actually 
creating a model file database. This unversal file is usually three times the memory file 
size of the .OP2 and model.POS file due to being in ASCII format. Model file database 
size (file.sfl and file.sf2) become approximately the same size as the universal file. 
Hence, a user should ensure plenty of hard drive memory is available during this 
translation phase, especially for large models (> 2000 elements). 
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IV. ONE-DIMENSIONAL ANALYSIS 



A. PLATE MODEL 



A one dimensional study was conducted utilizing the case defined by [Ref. 7] 
which involved modeling a floating plate on the surface of a bilinear fluid. The goal of 
this model was to accomplish the same results shown in [Ref. 7] by creating a finite 
element model and solving with the USA code. Follow-on study was performed by 
Newton [Ref. 20] in area of cavitation in support of developing a “shock shield,” a gas- 
filled cushion that was to surround the exterior of a submarine hull and reduce the 
magnitude of underwater shock loads transmitted through the hull. 

Fig. 10 and Table 1 describe the particulars of the Bleich-Sandler case. 



1 thin shell element, 4 nodes, 1 ” thick 
Symmetry on all four edges 
Dimensions: 1.5”x 1.5” 

E=30.0x 1 0 6 psi 

Poisson’s ratio=0.3 

Mass density=5 .32986 x lO^lb/in 3 

Table 1. Plate Characteristics 



Fluid Depth: 150.0 in 

100 volume elements, 8 nodes, 1.5”X1.5” 1.5” 
404 total fluid nodes 
Mass density= 9.3455 x 10' 5 lb s 2 /in 4 
Speed of sound= 5.712 x 10 4 in/sec 

Table 2. Fluid Characteristics 
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Figure 10. IDEAS Bleich Plate and Fluid 
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Peak pressure =103 psi 
Decay time= 9.958 E-4 sec 

Charge depth = 1 .0 E7 in, centered on the center of the plate 

Table 3. Incident Pressure Characteristics 



B. RESULTS 

Utilizing a time step size of 1.313 x 10" 5 sec, 1200 time steps were collected. 
With the CFA interface turned on with the USA, program Fig. 11 depicts the computed 
velocity and pressure at the plate. Fig. 12 depicts pressures at specific locations along 
the column of fluid. Because of the numerical computation scheme within the USA code, 
time zero noted on each pressure figure refers to the point in time that the incident 
pressure wave arrives at that particular fluid node and not the global time of plate motion. 
Total pressures at time zero for each fluid node is the sum of the atmospheric, 
hydrostatic, and incident pressure wave at that point of arrival from the initial charge 
location and applied exponential decay. Fig. 11 illustrates that the cavitation region does 
not touch the interface boundary between the fluid and the plate. Node 401 in Fig 12 
confirms the cavitation region is captured by the fluid FEM and does not extend beyond 
the DAA boundary. Fig. 13 depicts the plate velocity and pressure with CFA OFF. 

Calculated time periods of zero pressure or cavitation occurrence grow to a period 
of .008 sec half-way through the column and stay approximately the same to the end of 
the fluid column. In addition, the CFA captured the cascading effects of cavitation 
closure which occured at similar times following incident pressure wave arrivals at each 
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fluid node and closure pressures became larger at deeper depths than at the surfaces. This 
closure phenomenon begins at the point t=0108 sec at which the plate experiences an 
immediate velocity cutoff of the downward direction back to zero by the arrival of this 
“secondary shock” at the plate. Pressure at this particular time jumps to approximately 
14% of the peak pressure and creates an upward velocity of approximately 16% which 
brings the plate velocity back to zero. Multiplying time by 1004.22 to get it in terms of 
decay time units of the incident wave and normalizing velocity to 1/1000 of the speed of 
sound with a sign change due to the z-axis for verical displacement being defined 
downward into the fluid, calculated results match those (Fig. 14) with Bleich and Sandler 
as shown in Fig. 15 and 16. 

Hence, the finite element modeling approach with the USA/CFA code produced 
identical results as those shown by Bleich and Sandler through the method of 
characteristics [Ref. 7]. In addition to the pressure history profile of Fig. 15, USA 
postprocessing was able to display the cavitation closure anywhere in the fluid column. 
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Figure 11. Bleich Plate Velocity and Pressure (CFA ON) 
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Pressure(lbf/in A 2) Pressure(lbf/in A 2) Pressure(lbf/in A 2) 



DEPTH=3" (NODE 9) 




Figure 12. Bleich Fluid Pressure at Nodes 9, 200, and 401. (CFA ON) 
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Vertical Velocity 




Time(sec) 



Pressure 




Figure 13. Bleich Plate Velocity and Pressure (CFA OFF) 
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Normalized Vertical Velocity 




Figure 14. Normalized Vertical Velocity (— CFA ON, CFA OFF) 
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t , non-dimensional time 



12 - 




Figure 15. Time history of cavitated region [Ref. 7] 





Figure 16. Non-dimensional upward velocity of plate [Ref. 7[ 
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V. 



TWO DIMENSIONAL PROBLEM 



A. CROSS-SECTIONAL HULL 

The mid-section of the imported wireframe hull of the DDG model was utilized 
for this particular study, assuming structural properties remain linear throughout the 
process. Previous work conducted by Newton [Ref. 21] in the area of a two dimensional 
study involved a symmetric cylinder submerged underwater exposed to a plane wave 
incident shock front, utilizing earlier Fortran finite element programs for meshing, 
propagation of underwater shock, and transient response calculations. 

The goal in this particular model is to study responses, specifically peak 
velocities, with and without cavitation effects, of a surface ship model with a draft of 20 
feet. The time frame was limited to the first 30 msecs, in an attempt to eliminate the 
effects of the bubble pulse while comparing plane wave modeling and DAA2 spherical. 

Utilizing the USA code. Fig. 17 illustrates location of cross section in the global 
ship wireframe, 1 foot in depth, and Fig. 18 to Fig. 20 illustrate the meshed model. In Fig. 
18, red donotes the structure mesh while green denotes the fluid mesh. Fig. 20 depicts 
the outer DAA, hydrostatic, and structural interface elements mentioned previously. 
Tables 4-6 provide specific model characteristics. 



TOTAL NUMBER OF GRID POINTS 


452 


TOTAL NUMBER OF THIN SHELL ELEMENTS 


230 


THICKNESS OF ELEMENTS 


0.5 in 


TOTAL NUMBER OF WET ELEMENTS 


40 



Table 4. Structure Characteristics 
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MODULUS OF ELASTICITY (ksi) 


29600.0 


DENSITY (lbm/ft 3 ) 


489.02 


YIELD STRENGTH (ksi) 


51.00 


POISSON’S RATIO(v) 


.3 



Table 5. Element Properties 



NUMBER OF FLUID NODES 3034 

NUMBER OF FLUID VOLUME ELEMENTS 1 440 

NUMBER OF DAA FACE ELEMENTS 40 

NUMBER OF HYDROSTATIC PRESSURE FACE ELEMENTS 72 
NUMBER OF PRESCRIBED PRESSURE FACE ELEMENTS 0 
FACE ELEMENTS CONTACTING STRUCTURE 40 

MAXIMUM NUMBER OF FLUID VOLUME NODES 8 

MAXIMUM NUMBER OF FACE ELEMENT NODES 4 



Table 6. Fluid Characteristics 



Utilizing the bulk cavitation program to determine the extent of various bulk 
cavtiation zones of a TNT explosion, the smallest region of bulk cavitation takes place 
directly above the charge at the surface. Therefore, a charge location directly beneath the 
hull at 100 feet standoff (depth = 120 feet) would minimize effects of the bulk cavitation 
zone in the general vicinity of the hull model and provide a means for setting the lower 
edge of the DAA boundary just beyond the cavitating region at a depth of 50 feet, 30 feet 
under the keel. 



B. RESULTS 

Fig. 21 to Fig. 30 record the velocities at the locations in the structure and time 
varying pressures at locations in the fluid volume region for the CFA analysis noted on 
Fig. 19. Fig. 31 to Fig. 33 provide a comparison of keel velocities and accelerations. The 



44 



pressure profiles, like that seen in the one-dimensional analysis, reflect a respective time 
zero which starts at the point of incident pressure wave arrival at that particular node and 
not the initial time of structure motion. Structure motion is derived from the cumulative 
integration of all the pressure profiles in the fluid region superimposing each pressure 
profile at different spatial locations in time. 

The introduction of a fluid mesh beneath the structure and the turning on CFA 
captured the pressure varying fluid under the hull (Fig. 21 and 24) which significantly 
effected the velocity response of the hull. The cavitation effects on the hull during plane 
wave and DAA2 calculations produced almost identical results in the early time phase of 
hull velocity at the keel (node 42) which was the closest to the charge location. As 
expected in late time response, the plane wave solution begins to slightly vary from the 
DAA2 calculations. However, peak velocity response are the same in early time. 

Pressure profiles at the keel for both CFA models (Fig. 21 and 24), which are 
practically identical, indicate the first signs of cavitation closure pressure occurring at a 
much earlier time compared to the Bleich-Sandler case with an appreciably larger 
pressure peak, 41% of the initial pressure wave. Numerous other cavitation closure 
events occur at the keel with decreasing max pressures. This observed behavior was also 
an event reported in the DDG-53 USS John Paul Jones Live Fire and Testing Shock Post 
report [Ref. 22] which may have contributed significantly to large accelerations at 
particular electrical switchboards following the incident pressure wave. 

Keel velocities in the plane wave (CFA OFF) and DAA2 (CFA OFF) analysis, 
which places the DAA boundary at the ship’s wet hull surface, did not capture the 



45 



detailed high frequency response of the structure exposed to identical underwater 
explosion conditions (Fig. 27 and 29). Peak keel velocities were identical in magnitude 
and occurrence in time but greatly diverged past reaching the peak velocity of 162 in/sec. 
Pressures along the keel during both CFA OFF calculations were found to include only 
the large initial pressure wave arrival and small variations in the pressure later in time 
along the wet hull, while those in the DAA2 calculation being slightly larger. 

Internal deck velocities, however, revealed velocities greater during the plane 
wave modeling with a surrounding fluid mesh and CFA ON (Fig. 22). The pressure and 
corresponding variations along the wet surface of the hull caused internal velocities to be 
greater than those seen in the DAA2 CFA ON solutions (Fig. 25), primarily due to the 
less accurate calculations late in time for the plane wave equations and higher velocity 
response calculations. Both CFA OFF runs (Fig. 28 and 30) reveal small variations in 
velocity, the greatest occuring at the lowest deck. These small velocites could be 
attributed to the simplified cross sectional hull model which does not include any internal 
longitudinal or transverse bulkheads. If these particular bulkheads were modeled between 
decks, the transit velocity responses of the keel could have been transmitted to the upper 
decks instead of the bending the outer hull only. 

Fig. 31 depicts the overlay of velocities of all 4 models (plane wave with and 
without CFA, DAA2 with and without CFA) at the keel, illustrating the larger velocities 
of the CFA OFF calculations for both plane wave and DAA2 modeling. Plane wave 
modeling of the underwater explosion incident pressure wave without a fluid mesh 
provides for a conservative approach in calculating the peak velocity of the structure at 
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the fluid structure interface. A shift to a higher velocity at an earlier time was more 
pronounced during this particular CFA study in comparison to the Bleich-Sandler case 
with a 22 % increase in peak velocity. Appreciable differences in acceleration were noted 
between CFA ON and CFA OFF calculations (Fig. 32 and 33). CFA OFF results 
reflected initial accelerations 3.5 times larger than CFA ON results but with no following 
variations. 

Peak velocities in each case are provide in Table 7. 



LOCATION 


Node 


Plane wave 
(CFA ON) 


DAA2 
(CFA ON) 


Plane wave 
(CFA OFF) 


DAA2 
(CFA OFF) 


Keel 


42 


124 


124 


162 


162 


Inner Bottom 


3209 


Increasing 


Increasing 


Increasing 


4.3 


Deck 3 


3234 


Increasing 


Increasing 


.094 


.13 


01 Level 


3162 


Increasing 


Increasing 


2.7E-4 


Increasing 



Table 7. Maximum Velocities (in/sec) 



The responses labeled “increasing” in Table 7 reflect increasing values at the end of the 
time period studied. 

Fluid pressures at nodes 3021, 1600, and 2994 during both plane wave and DAA2 
CFA ON calculations (Fig.23 and 26) show similar relative patterns of growing numbers 
and durations of cavitation periods as depth decreases. The magnitude of peak pressures 
following initial cavitation in the DAA2 results tended to be 15-20% larger than the plane 
wave peak pressures. Node 3021 pressures during both CFA ON calculations at t=.0025 
secs captured the largest peak pressure, followed by Node 2994 at t=.0047 secs. 
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Hence, based on the results, the best conservative approach to modeling a surface 
ship exposed to an underwater explosion would be to apply the plane wave incident 
pressure wave in a fluid mesh that surrounds the entire wet surface of the ship. 
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Figure 18. Structure and Fluid FEM 
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Figure 19. Cross Section Dimensions 




Figure 20. DAA, Hydrostatic, and Wet Structure Boundaries 
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Figure 21. Velocity and Pressures at the Keel (Node 42)-Plane wave, CFA ON 
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x 10 



Node 3162 






Figure 22. Internal Deck Velocities-Node 3162(top), Node 3234 (middle), Node 3209 

(bottom). Plane wave, CFA ON 
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Pressure(lbf/in A 2) Pressure(lbf/in A 2) 



Node 3021 (DAA Boundary) 




Figure 23. Fluid Pressures-Node 3021 (top), Node 1600 (middle), Node 2994 

(bottom). Plane wave, CFA ON 
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Vertical Velocity 




Pressure 




Figure 24. Velocity and Pressure at the Keel (Node 42)- DAA2, CFA ON 
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Figure 25. Internal Deck Velocities- Node 3162 (top), Node 3234 (middle), Node 3209 

(bottom). DAA, CFA ON 
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Pressure(lbf/in A 2) Pressure(lbf/in A 2) 



Node 3021 (DAA Boundary) 




Node 1600 (Mid-Fluid) 




Figure 26. Fluid Pressures-Node 3021 (top), Node 1600 (middle), Node 2994 

(bottom). DAA2, CFA ON 
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Pressure 




Time(sec) 

Figure 27. Velocity and Pressure at the Keel (Node 42)-Plane wave, CFA OFF 
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Figure 28. Internal Velocities-Node 3162 (top), Node 3234 (middle), Node 3209 

(bottom). Plane wave, CFA OFF 
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Figure 29. Velocity and Pressure at the Keel (Node 42)-DAA2, CFA OFF 
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Figure 30. Internal Velocities-Node 3162 (top). Node 3234 (middle), Node 3209 

(bottom). DAA2, CFA OFF 
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Figure 31. Comparison of Keel Velocities. 
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KEEL VELOCITY (NODE 42) 
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Planewave CFA OFF 




Figure 32. Comparison of Keel Accelerations-PIane wave 
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Figure 33. Comparison of Keel Accelerations- DAA2 
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VI. THREE DIMENSIONAL MODEL 



The ideal three dimensional ship analysis would include: 

1 . A detailed ship model which included specific hull and equipment 
characteristics (locations, weights, etc.) 

2. A fluid mesh model of the surrounding ocean that captures the cavitation region 
in three dimensions, as seen in the one and two dimensional analysis. 

3. A computer hardware that can handle the detailed processing and large solution 
database. 

Due to the complexity of a global ship model and the time integration required, a 
simplified “notional” ship design was utilized which was part of an earlier model of the 
Arleigh Burke Destroyer. Based on the results obtained with the two-dimensional model 
and the requirements for conditions 2 and 3 stated above, a plane wave/CFA OFF 
analysis was utilized for the three-dimensional ship model which sets the DAA boundary 
at the wet surface of the hull. Table. 8 provides general characteristics of the DDG model 
[Ref. 23]. 



Length Between Perpendiculars (LBP) 


466.0 ft 


Length Overall (LOA) 


492.1 ft 


Beam, Weather Deck 


66.5 ft 


Draft to Keel 


20.7 ft 


Lightship weight 


6327 LTON 


Full Load displacement 


8000.3 LTON 



Table 8. DDG-51 ASSET Design Summary [Ref. 23] 



Fig. 34 and 35 depict solid modeling of the imported DDG model and specific 
charge locations for the two models created. Fig. 36 illustrates locations of the following 
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major areas of interest were selected for their general purpose of holding vital pieces of 
equipment such as the AEGIS weapon system components and main propulsion units: 
node closest to the charge, bow, stem, Main Engineering Rooms (MER 1 and 2), Combat 
Systems Equipment Rooms (CSER 1 , 2, and 3), and the general location of th SPY radar 
room in the superstructure. 
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Figure 34. Solid Modeling of ASSET DDG-51. 
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Figure 35. Charge Locations Below the Keel (1) and Offset (2). 
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Figure 36. General Space Locations. 
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A. 



1962 FINITE ELEMENT MODEL 



1. Modeling 

This initial model utilized a general plate thickness of 0.4 inches of similar 
material qualities used in the two dimensional model throughout the entire ship. Inner 
bottom stiffeners and heavy equipment locations were not constructed in the original 
IGES model out of ASSET but were later created in IDEAS for the 5682 element model 
study. Fig. 35 illustrates the 1962 element model which includes 348 wet elements. In 
order to capture peak velocities, damping coefficients were not applied to the model. 

Explosive charge conditions were the same as those utilized for the two 
dimensional analysis: 



CHARGE TYPE 


TNT 


WEIGHT 


100 LBS 


LOCATION (relative to ship) 


MIDSHIPS, CENTERLINE, 100 ft beneath the keel 


DEPTH 


120 feet 



Table 9. Charge characteristics 



2. Results 

The whipping motion of this particular model exhibited the expected basic motion 
of a simple free-free beam excited by a forcing function in the middle of the model. The 
highest velocities and accelerations occurred at the node (Node 42) closest to the charge 
(Fig. 37 and 38) due to an unsmeared and unstiffened keel and outer hull. The 
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unsmeared and unstiffened keel also experienced high frequency oscillations in 
comparison to other areas of the model. The bow and stem experienced the second 
largest maximum vertical velocities and accelerations (Fig. 39 to 42) which occurred 
approximately 45 msec and 32 msec, respectively, after the peak keel velocity. 

Peak velocities at the bow, stem, and closest point to the charge at the keel indicate 
the greatest magnitude in velocities to be in the vertical direction, followed by fore and 
aft motion, and the athwartships direction to be the lowest. However, internal deck 
locations exhibited the greatest peak velocities in the fore and aft direction, followed by 
athwartship, and vertical velocities to be the lowest because of their general locations 
close to the incident shock wave arriving midships in the model. Both MER spaces (Fig. 
43 to 46) experienced the largest internal accelerations, specifically in the fore and aft 
direction followed by CSER 3 (Appendix C) with almost equal velocites and 
accelerations in the fore and aft and athwartship directions. MER 2 experienced the 
largest fore and aft motions in comparison to the other internal spaces studied, largely due 
to its close location to the initial shock wave. Amongst the combat systems spaces, the 
SPY radar room (Appendix C) had the highest fore and velocity (9 in/sec) but with 
similar accelerations. The asymmetry of the aft super structure arrangement places a 
majority of the structure to port. This asymmetric condition in the DDG-51 model 
contributed to the athwartship velocities and accelerations in all areas of the ship. 

Modeling the ship globally without a fluid volume mesh and exposing it to a 
plane wave underwater explosion captured the high frequency vibration or “ringing” in 
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the hull which was not observed in the plane wave/CFA OFF two dimensional cross- 
sectional hull model. Table 10 summarizes maximum velocities and accelerations: 



Location 


Node 


Fore and Aft 


Athwartships 


Vertical 






V(in/sec) 


A(G) 


V(in/sec) 


A(G) 


V(in/sec) 


A(G) 


Keel 


42 


10.2 


47.1 


3.3 


20.8 


176 


2210 


Bow 


842 


13.1 


9.2 


9.1 


4.3 


27.2 


10.1 


Stem 


2157 


13.4 


10.1 


3.8 


1.33 


20.1 


12.2 


MER1 


1933 


7.2 


12.1 


2.6 


2.2 


.7 


.058 


MER2 


1895 


12.3 


9.92 


4.6 


3.8 


4.2 


.13 


CSER1 


1649 


5.5 


3.6 


2.5 


1.3 


.17 


.019 


CSER2 


1748 


6.5 


5.1 


1.33 


1.1 


.16 


.025 


CSER3 


1990 


6.1 


7.1 


6.8 


5.03 


.135 


.08 


SPY Radar 
Room 


2534 


9 


4.8 


2.6 


1.8 


.055 


.0088 



Table 10. Maximum Absolute Velocity and Acceleration (1962 Elements) 
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Velocity (in/sec) . Veloci >'' < in/sec > 



Fore and Aft Direction 




Time (sec) 
Athwanships 




Time(sec) 

Vertical 




Time (sec) 



Figure 37. Keel Velocity (Node 43) 
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Acceleration (G) 



Fore and Aft Direction 





Time(sec) 

Figure 38. Keel Acceleration (Node 43) 
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Velocity (in/sec) Velocity (in/sec) Velocity (in/sec) 



Fore and Aft Direction 






Time(sec) 



Figure 39. Bow Velocities (Node 842) 
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Acceleration (Q) Acceleration (G) , Acceleration (G) 



Fore and Aft Direction 






Figure 40. Bow Accelerations (Node 842) 
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Velocity (in/sec) Ve|ocjty (in/sec) Velocity (in/sec) 



Fore and Aft Direction 




Timetsec) 

Athwartships 




Figure 41. Stern Velocities (Node 2157) 
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Acceleration (G) Acceleration (G) , Acceleration (G) 



Fore and Aft Direction 





Time(sec) 

Vertical 




Figure 42. Stern Accelerations (Node 2157) 
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Fore and Aft Direction 





Time(sec) 

Vertical 




Time(sec) 

Figure 43. MER 1 Velocities (Node 1933) 
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Fore and Aft Direction 





Figure 44. MER 1 Accelerations (Node 1933) 
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Velocity (in/sec) Velocity (in/sec) . Velocity (in/sec) 



Fore and Aft Direction 





Figure 45. MER 2 Velocities (Node 1895) 
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Fore and Aft Direction 





Time(sec) 



Figure 46. MER 2 Accelerations (Node 1895) 
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B. 5682 FINITE ELEMENT MODEL 



1. Modeling 

In order to provide a detailed mesh for the wet surface of the hull and improve 
equipment locations, a finer mesh was created of the same model which included inner 
bottom stiffeners and 955 wet elements, Fig. 47. Utilizing the available ASSET data on 
main engineering equipment [Ref. 23], early Flight II/IIA study drawings [Ref. 24], and 
DDG combat system equipment data [Ref. 25], large groups of heavy equipment were 
lump massed in their respective areas. 

The detailed requirement for longitudinal and transverse stiffners was simplified 
by the smearing technique utilized in ASSET (APPENDIX D). Although previous work 
by Cunningham [Ref. 26] revealed the limitations in local and global smearing, the 
ASSET technique sizes stiffeners to provide a T-section with reasonable proportions that 
meet strength and geometric constraints. The stiffener cross-sectional area is smeared into 
the plate cross sectional area, creating a smeared plate thick than the standard plate by a 
factor of (1 + smear ratio(S)). Based on the data for the notional DDG design which 
included intitial plate thickness and stiffener smearing ratios, the following plate 
thicknesss were used in this particular model: 
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Structure Location 


Smeared Plate Thickness (in) 


BOTTOM SHELL 


.57 


INNER BOTTOM 


.42 


DECK 1 


.51 


DECK 2 


.54 


DECK 3 


.49 


DECK 4 


.48 


MAIN DECK 


.55 


TRANSVERSE BULKHEADS 


.33 


OUTER HULL 


.47 


SUPERSTRUCTURE 


.33 



Table 11. Smeared plate thicknesses (inches) 



In addition to Appendix E which denotes general weights and locations of 
combat system equipment that were lump massed, the following specific engineering 
room equipment was used with respect to their locations derived from [Ref. 24]. 



GAS TURBINES (4) 


84.7 


GAS GENERATORS (3) 


94.2 


REDUCTION GEARS (2) 


109.5 



Table 12. Engineering Equipment Total Weights (LTON)[Ref . 23] 



The net weight of the entire model was 3598 Lton, 43% lighter than the original 
“light weight” DDG ASSET model. 

Charge location was placed in an offset position utilizing only 120 lbs TNT at a 
standoff distance of 107 feet and a depth of 50 feet . 
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2 . 



Results 



Due to the increased complexity of the meshed model, post processing of this 
particular model . became quite difficult during the MSC/NASTRAN restart which 
prevented further post-processing into IDEAS. Fig. 48 to 51 depict the greatest velocities 
and accelerations occuring midships of this model, specifically at the point along the hull 
closest to the charge (Node 393) and at the node along the keel with the same axial 
distance from the bow (Node 397). These two locations coincide with those same points 
on the hull noted previously in Fig. 35 with charge 1 for node 397 and charge 2 for node 
393. 

Bow and stem reactions (Fig. 52 to 55) were noted higher than those seen in the 
1962 element model due to the majority of concentrated masses being located in central 
locations of the ship (i.e. aft of the forward VLS and forward of CSER3) as well as the 
spaces of interest (Fig. 36) and the lack of masses in the general area of the bow and 
stem. Peak bow and stem vertical reactions occurred 40 msec and 30 msec, respectively, 
later relative to midships reactions (Node 393 and 397) with high frequency oscillations. 
MER 1 (Fig. 56 and 57) exhibited lower velocities and accelerations relative to MER2 
(Fig. 58 and 59) which can be contributed by the asymmetry between the two 
engineering rooms. The gas turbine modules and reduction gear masses were located on 
the port side of the engine room of MER 2 which was the same side as the charge 
location and, as seen in the 1962 element model, MER 2’s location was closer to the 
incident shock wave arrival point along the hull. Bow,stem, and both MER spaces tended 
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to have reactions in the following order of highest to lowest magnitude: vertical, fore/aft, 
and athwartships. Combat system equipment rooms (CSER 1,2 and 3) and the SPY radar 
room (Appendix F) showed very little motion in the fore/aft and athwartships direction 
with practically zero vertical motion. 

For this particular tactical scenario and model, the only shipboard internal space 
which experienced the highest acceleration was MER 2 with large magnitudes of 
oscillations in the fore and aft direction as well as athwartship (Fig. 58 and 59). All other 
internal spaces encountered low accelerations which did not seriously introduce severe 
shock loadings upon their respective equipment. The effect of smearing and lump 
massing groups of large pieces of internal equipment created dramatically different 
results compared to the 1962 unsmeared, unlumped element model. In addition to these 
structural differences, the low charge weight and offset distance of the underwater 
explosion did not produce appreciable effects on internal ship behavior. Table 12 
provides a summary of the maximum velocities and acclerations in each area: 
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Location 


Node 


Fore and Aft 


Athwartships 


Vertical 






V(in/s) 


A(G) 


V(in/s) 


A(G) 


V(in/s) 


A(G) 


Standoff 

Point 


393 


17.1 


64.1 


53.2 


172 


58.9 


210 


Keel 


397 


10 


52 


30.2 


110 


28.8 


106 


Bow 


1554 


13.2 


42.3 


12.1 


44 


21.2 


58 


Stem 


2281 


13.8 


23.1 


10.8 


30.1 


12.2 


35.6 


MER 1 


4448 


.33 


.18 


2.1 


.67 


2.23 


1.08 


MER2 


4543 


8.8 


42.4 


9.2 


44.5 


39.3 


133 


CSER1 


3567 


.013 


1.9E-3 


.012 


1.3E-3 


6.4E-4 


5.1E-5 


CSER2 


3669 


.19 


.023 


.78 


.07 


7.7E-4 


1.48E-4 


CSER3 


3541 


.15 


.19 


1.4 


.23 


7E-3 


1.5E-3 


SPY Radar 
Room 


6955 


.59 


.071 


.15 


.019 


3.2E-4 


4.3E-3 



Table 13. Maximum A 



asolute Velocities and Accelerations (5698 element model) 
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Velocity (in/sec) Velocity (in/sec) Velocity (in/sec) 



20 



Fore and Aft Direction 




Timefsec) 

Athwartships 




Time(sec) 



Figure 48. Hull Standoff Point Velocities (Node 393) 
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Acceleration (G) 



Fore and Aft Direction 





Time(sec) 

Vertical 




Figure 49. Hull Standoff Point Accelerations (Node 393) 
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Velocity (in/sec) Velocity (in/sec) Velocity (in/sec) 



Fore and Aft Direction 





Figure 50. Keel Velocities (Node 397) 
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Fore and Aft Direction 





Figure 51. Keel Accelerations (Node 397) 
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Velocity (in/sec) Velocity (in/sec) Velocity (in/sec) 



20 



Fore and Aft Direction 





Figure 52. Bow Velocities (Node 1554) 
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Fore and Aft Direction 




Time(sec) 



Figure 53. Bow Accelerations (Node 1554) 
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Velocity (in/sec) Velocity (in/sec) Velocity (in/sec) 



Fore and Aft Direction 





Time(sec) 

Vertical 




Figure 54. Stern Velocities (Node 2281) 
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Acceleration (G) Acceleration (G) Acceleration (G) 



Fore and Aft Direction 





Time (sec) 
Vertical 




Time(sec) 



Figure 55. Stern Accelerations (Node 2281) 
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Velocity (in/sec) Velocity (in/sec) 



Fore and Aft Direction 




Figure 56. MER 1 Velocities (Node 4448) 
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Acceleration (G) Acceleration (G) Acceleration (G) 



Fore and Aft Direction 




Time(sec) 

Vertical 




Figure 57. MER 1 Accelerations (Node 4448) 
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Velocity (in/sec) Velocity (in/sec) Velocity (in/sec) 



Fore and Aft Direction 





Figure 58. MER 2 Velocities (Node 4543) 
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Fore and Aft Direction 





Figure 59. MER 2 Accelerations (Node 4543) 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



The introduction of a fluid volume model with a structure and utilizing the USA 
/CFA interface provides a means to capture the effects of cavitation upon a surface vessel 
in the early time frames of an underwater explosion. The effects of cavitation closure as 
seen in the two dimensional model can be quite significant and inflict another impulsive 
force upon the ship as reported in the DDG 53 shock trial reports [Ref. 22] in addition to 
the incident wave, bubble pulsations, bottom reflection, or bottom refraction. This 
increases the complicated task of predicting surface ship behavior when exposed to an 
underwater explosion. CFA OFF modeling may provide a conservative approach in 
predicting ship response to an underwater explosion in early time. Finite element 
modeling, time integration calculations, and tactical placement of an explosive device 
must be carefully done to accurately model surface ship behavior. Plane wave modeling 
with a surface ship model and a fluid mesh may provide more accurate overall ship and 
internal equipment responses. 

With improving systems and capacity, the testing of new ship models with 
specific interior designs from the ship designer with the underwater shock analyst 
utilizing the USA/CFA code can prove to be a formidable method in predicting forcing 
functions on internal equipment of a surface vessel. In view of the changes of the DDG 
Flight IIA designed compared to the present DDG Flight II, this particular method may 
provide some cost benefits in future ship shock trials. 
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Recommend continuing studies in the following areas: 

1 . Improve on the two and three dimensional model with internal masses or fluids. 

2. Obtain actual DDG-51 Flight IIA DDG designs from NAVSEA which can be 
imported into an FEM program such as IDEAS or MSC/PATRAN which can be 
tested with the USA code for both linear and non-linear responses. 

3. Introduce other surface ship model types such as the slice hull. 

4. Model the underwater tactical scenario with significant shallow water bottoms 
in order to study “littoral” responses of a surface ship close to the shore. 
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APPENDIX A. BULK CAVITATION PROGRAM 



The following program calculates the bulk cavitation zone for various weights and 
depths of TNT. The program can be easily modified for other explosive material. 
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%The following program will compute the bulk cavitation zone of an underwater 
%explosion of TNT explosive with depth and weight of charge as inputs. 

% 

%TNT characteristics 

Al=l . 18 ; A2=- . 185 ;K1=22505 ; K2= . 058 ; 



%D= input (' Enter depth of your charge ') 
%W=input (' Enter the weight of TNT explosive ') 
Co=4967 ; 

Pad 4 . 7 ; 

x= (0 : 5 : 1200 J ' ; 
n=length(x) ; 
yu=diag ( zeros (n) ) ; 
yl=diag ( zeros (n) ) ; 
for j=l:n 

y(l) =0; 

for i=l:500 

rl = sqrt ( ( (D-y (i) ) ~2) + (x { j ) ~2) ) ; 
r2=sqr t ( ( (D«-y(i))*2)«-(x(j)"2)); 
theta=K2*W~ (1/3) * ( (W~ (1/3) /rl) ~A2 ) ; 
dr = (r2-rl) / (Co* theta/ 1000 ) ; 

Pmax=Kl * ( (r(l/3)/rl) A Al) ; 

Pl=Pmax*exp { -dr ) ; 

Pi =Pmax*exp ( -dr ) ; 

P2=K1* ( (VT (1/3) /r2)*Al) ; 

Ph=64 *y ( i ) /144; 



%speed of sound in water (ft/sec) 
%atmospheric pressure 



%real charge distance 
%image charge distance 
%decay constant (msec) 

%decay rate 
% Incident shock wave 

%Inc shock wave 0 raref arriv 
%Inc shock wave Oraref arr 
%rarefrac at arrival time 
%hydrostatic pressure 



F=Pl+Pa+Ph-P2 ; 

if ( F< = 0 ) 

break 

else 

y ( i ♦ 1 ) =y(i) + . 1; 

end 

end 

yu ( j ) =y ( i ) ; 
y(l) =0; 

for s = 1 : 500 

r l = sqrt ( { (D-y (s) ) ~2 ) + (x ( j ) ~2 ) 
r2=sqrt ( ( ( D+y (s) )~2)+(x(j)~2)) ; 
theta=K2*W~ (1/3) * ( (VT (1/3 ) /rl ) ~A2 ) ; 
Pmax=Kl * ( (W~ (1/3) /rl) ~A1 ) ; 

Pi = Pmax * exp ( -dr ) ; 

P2 =K1 * ( (W~ (1/3) / r2 ) Al ) ; 

Ph=64 *y ( s ) / 1 4 4 ; 

dr= ( r2 -r 1 ) / (Co* theta/ 1000 ) ; 

gl= ( D+y ( s ) ) / r2 ; 

g2=r2- (2*D*gl) ; 

g3=(g2/rl) * ( ( (A2*r2) /rl) -A2-1) ; 
g4 = Pi / (Co* theta/1000) ; 
g5= { Al * Pi ) / (rl~2) ; 
g6=Pl+Pa+Ph; 

G=-(g4*(l+g3) )-(g5*g2) +(64*gl/144)+ 
if (G>=0 ) 
break 
else 

y(s+l)=y(s)+. 1 ; 

end 

end 

y 1 ( j ) =y ( s ) ; 

if (yl ( j ) <=yu ( j ) ) 
z = j ; 
break 

else 

z = n ; 

end 



%F(x,y) equation for upper boundary 



%real charge distance 
%image charge distance 
%decay constant (msec) 

% Incident shock wave 
%Inc shock wave @ raref arriv 
%rarefrac at arrival time 
%hydrostatic pressure 
%decay rate 



(Al*g6/r2); %Combined G(x,y) lowcav bound 



end 

yu = ( - 1 ) * yu ( 1 : z ) ; 
y 1 = ( - 1 ) * y 1 ( 1 : z ) ; 

X f = x ( 1 : z ) ; 
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APPENDIX B. CROSS-SECTIONAL HULL INPUTS 



The following are sample inputs for the two-dimensional model only. 



assign outputs * uti ' , umc=u , delete $ 
assign output2= ' IDEAS1 . op2 ' , unit=12 , unknown 
ID USA, CROSS HULL 3 
SOL 109 $ MSC 

DIAG 8,13 $ 

TIME 99 

$ 

$ START SDRC'S I-DEAS OPTION 

compile phaseO souin=mscsou nolist noref $ 

$ add CALL to get output2 geometry datablocks early in solution, 
alter 297 $ before call seprep2 
type parm, , i, y, ounit2=12 $ 
if (post < 0) then $ 

OUNITl=OUNIT2 

CALL OUT2GEOM CSTM, GPL, GPDT, EPT, MPT, GEOM2 , GEOM3 , GEOM4 , 
GE0M1,DIT, DYNAMICS, EQEXIN, BGPDT , CASECC// 
OUNITl/O $ 

endif $ 

$ END SDRC'S I-DEAS OPTION 

$ 

compile subdmap=selg , souin=mscsou, noref , nol ist $ 

ALTER 14 $ AFTER SSG1 

type parm, , i , n, storit = 0 

call dbstore pjx, ,,, //seid/0/' DBALL' /s, storit $ 

OUTPUT2 BGPDTs , EQEXINs , , ,///ll//25000 $ 

$ 

compile subdmap=sekr, souin=mscsou, noref , nolist $ 

ALTER 49 $ AFTER gpsp 

type parm, , i , n, foundit = 0 $ 

call dbfetch /pjx, ,,, /seid/0/0/0/s, foundit $ 

OUTPUT2 uset , p j x, cstms , , ///11//25000 $ 

$ 

compile subdmap=DTRANRS, souin=mscsou, noref , nolist $ 
alter 1 
type db ECTS 

ALTER 6 $ before TRD1 

OUTPUT2 ECTs, BDD , KDD , MDD, ///11//25000 $ 

$add EOF marker 
OUTPUT2 //- 9/12/ $ 

$ 

EXIT $ 

CEND 

TITLE=CROSS HULL 

LOADSET=10 

TSTEP=25 

SPC = 1 

DISP (plot ) =ALL 

velocity (plot ) =all 

stress (plot) =all 



BEGIN BULK 

PARAM AUTOS PC YES 

PARAM POST - 2 



$ TSTEP, SID, Nl, DTI , NOl 
TSTEP, 25,40,0.025,1 
$ 

$ LSEQ , S I D , DAREA , L I D 
LSEQ, 10, 10,1 

GRID 1 0 -336.0 

GRID 2 0 -336.0 



240.0 12.0 0 

240.0 0.0 0 
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CROSS SECTIONAL HULL3B PROBLEM 

1 , 0 , 0 , 1 , 1 , 0, 0 



USAF 

GAL 

F F F F 
F F F 



AUGMAT INPUT FOR CROSS HULL3 PROBLEM 
GAL HULL 3 . FLU HULL 3 . GEO HULL3 . PRE 



$ STRNAM FLUNAM GEONAM PRENAM 



F F F F 


$ 


FRWTGE 


FRWTST 


FRWTFL 


LUMPFM 


F F F T 


$ 


FLUSKY 


DAAFRM 


SYMCON 


DOFTAB 


T T F F 


$ 


PRTGMT 


PRTTRN 


PRTSTF 


PRTAUG 


F F F F 


$ 


MODTRN 


STRLCL 


INTWAT 


CFADYN 


NASTRAN-MSC 


$ 


MAINKY 








0 


$ 


NTYPDA 








452 1356 6 3 


$ 


NSTR 


NSFR 


NFRE 


NFTR 


1 


$ 


NSETLC 








0 1 80 1 


$ 


NDICOS 


J START 


JSTOP 


JINC 



FLUID MASS RUN FOR CROSS SECTION HULL3B 
HULL 3 . FLU HULL3.GEO GAL HULL3 . DAA 

T T T T 
T T F F 
F F T T 
F F F F 
F T F F 
F F F F 
F T F F 
F F 



FLUID. VOL FLUID. VOL 
NASTRAN-MSC 
0 452 0 80 
0 0 40 
0 3034 

9.3455-5 57120. 

10 

240 . 0 . 1 . 0 . 

14 . 7 386 . 4 
0 
0 
0 



$ FLUNAM GEONAM GRDNAM DAANAM 
$ PRTGMT PRTTRN PRTAMF CAL CAM 
$ EIGMAF TWODIM HAFMOD QUAMOD 
NASTAM STOMAS STOINV 
FRWT-GE FRWTGR FRESUR 
STOGMT ROTGEO ROTQUA 
STRMAS SPHERE ROTSYM 
CAVFLU FRWTFV INTCAV 
MASREF 
FACNAM 



PCHCDS 
FRWTFL 
RENUMB 
PRTCOE 
OCTMOD 
BOTREF 
CAVNAM 
MAINKY 
$ NSTRC NSTRF 
$ NBRA 
$ NFVNC 
$ RHO 
$ NVEC 
$ DEPTH 
$ PATM 
$ NUMCON 
$ NSRADI 
$ NSORDR 



NCYL 

NFVNF 

CEE 

CXFS 

GRAVAC 



NGEN 

NCAV 



CYFS 



NGENF 



C2FS 
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USA-NASTRAN CFA TIMINT RUN FOR CROSS HULL3 PROBLEM 



HULL3.PRE HULL 3 


. POS 


$ PRENAM POSNAM 




HULL3 . RES 




$ RESNAM WRTNAM 




F F F 


$ 


REFSEC 


FLUMEM 


STRVEL 




1 0 


$ 


NT I NT 


NCHGAL 






0. 3.E-5 


$ 


STRTIM 


DELTIM 






3.E-2 0. 0. 


$ 


FINTM 


ADAMP 


BDAMP 




T F F F 


$ 


EXPWAV 


SPLINE 


VARLIN 


PACKET 


F T F F 


$ 


HYPERB 


EXPLOS 


DOUBDC 


VELINP 


F F F 


$ 


BUB PUL 


SHKBUB 


BOTREF 




1 


$ 


NCHARG 








0. 14.7 


$ 


HYDPRE 


PATM 






0. -1.44E3 0. 


$ 


XC 


YC 


ZC 




0. -351.8 0. 




$ sx 


SY 


sz 




201 


$ 


JPHIST 








1. 0. 


$ 


PNORM 


DETIM 






3.E-5 


$ DTHIST 








2 


$ 


CHGTYP 








100. 100. 120. 




$ WEIGHT SLANT CHGDEP 


1 150 0 


$ 


NSAVER 


NRESET 


NSODFL 




0 0 0 0 


$ 


LOCBEG 


LOCRES 


LOCWRT 


NSTART 


T F F 


$ 


FORWRT 


STBDA2 


ASCWRT 




T F F 


$ 


NOAMBI 


PRTVOL 


PRTINI 




.5 0. 


$ 


FVBETO 


FVBET1 






Ill 


$ 


ICAVSW 


IORDF 


IORDV 




124 1440 3034 


$ 


NFVWAV 


NFVELM 


NFVNOD 




T 


$ 


DISPLA 








0 0 


$ 


NPREVT 


NPREVF 






T T 


$ 


LISTRE 


PRTPLT 






2 0 0 


$ NWETHS NDRYHS NUMSET 




42 2 


$ 


NODOUT 


NFROUT 






3209 2 


$ 


NODOUT 


NFROUT 






2 0 0 


$ 


NWETHS 


NDRYHS 


NUMSET 




42 2 


$ 


NODOUT 


NFROUT 






3209 2 




$ NODOUT NFROUT 




4 0 


$ 


NPREHS 


NUMSET 






42 0 


$ 


NEQHPR 


IPRS 






1600 1 


$ 


NEQHPR 


IPRS 






2994 0 


$ 


NEQHPR 


IPRS 






3021 0 


$ 


NEQHPR 


IPRS 






F 


$ 


SCALEF 









hull3B MODEL 
HULL3 . POS 
1001 , 40,0,0,1 
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assign inputt2= ' UT1 ' , unit=ll $ 
assign master= ' hull3B . MASTER' $ 
assign output2=' IDEAS 2 . op2' , unit =12 , unknown 
RESTART version=l / keep $ 

ID USA, HULL3B 
SOL 109 $ MSC 

DIAG 8,13 $ 

TIME 120 

$ 

compile subdmap=DTRANRS, souin=mscsou, noref $ 
ALTER 7,7 $ eliminate TRD1 
delete /PNLD, , , , / $ 

INPUTT 2 /UHV, , , , //ll $ 

EQUIVX UHV/UDVT/-1 $ 

CEND 

TITLE= THESIS HULL 3 B MODEL 

LOADSET=10 

TSTEP=4 6 

DISP (plot ) =ALL 

STRESS (plot) =ALL 

VELOCITY (plot) =ALL 

ACCELERATION (plot) =A11 

BEGIN BULK 

TSTEP ,46, 1001, 3 . Oe-5, 1 
ENDDATA 



108 



APPENDIX C. COMBAT SYSTEMS AREA RESPONSES-1962 ELEMENT 

MODEL 



The following plots represent the responses for CSER 1,2, 3 and the SPY radar 
room for the 1962 element model. 



/* 
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Figure 60. CSER 1 Velocities (Node 1649) 
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Figure 61. CSER 1 Accelerations (Node 1649) 
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Fore and Aft Direction 




Time(sec) 

Figure 62. CSER 2 Velocities (Node 1748) 
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Figure 63. CSER 2 Accelerations (Node 1748) 
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Figure 64. CSER 3 Velocities (Node 1990) 
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10 



Fore and Aft Direction 
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Figure 65. CSER 3 Accelerations (Node 1990) 
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Fore and Aft Direction 




Figure 66. SPY Radar Room Velocities (Node 2534) 

c5 
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Fore and Aft Direction 





Figure 67. SPY Radar Room Accelerations (Node 2534) 
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APPENDIX D. SMEAR PLATE DATA 

The following illustrates the smeared plate area method and ratios that were 



extracted from [Ref. 23]. 
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Figure 68. Smeared Plate Area Method [Ref. 14] 
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1 


11.00 


.66 


2 
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. 72 


3 


8.01 


. 71 


4 


10.01 


1.06 
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IN 

•51 28.56 

•69 19.71 

•68 19.78 

•95 36.86 



SMEAR 

FLANGE WT/FT RATIO 
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2.96 37.37 .10 

3.86 27.61 .22 

3.84 27.22 .22 

6.64 33.99 .22 
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APPENDIX E. COMBAT SYSTEM EQUIPMENT WEIGHTS 



The following table comprises the general combat system equipment utilized in 
the distribution of lump massing in the 5698 element model whose weights derived from 
[Ref. 25] and placed in locations utilizing [Ref. 24]. 



Space 


Location 


Weight (lb) 


Communication Center 


2-126-1-C 


1200 


Radio Trans. Room 


2-158-1-C 


5884 


Array Room No. 1 


03-128-2-Q 


14692 


Array Room No. 2 


03-128-2-Q 


14692 


Array Room No. 3 


03-155-1-C 


14692 


Array Room No. 4 


03-1 55-2 -C 


14692 


Radar Room No. 1 


03-128-0-C 


27720 


Radar Room No. 2 


03-142-0-C 


27720 


Radar Room No. 3 


01 -274-1 -C 


12035 


Power Supp/Conv Room 


3-126-0-Q 


11850 


Power Conversion Room 


3-319-0-Q 


11850 


Sonar Equip. Room No. 1 


1-18-0-Q 


11545 


Sonar Equip. Room No. 2 


2-18-0-Q 


11545 


Sonar Equip. Room No. 3 


3-18-0-Q 


11545 


Sonar Dome 


External to bow 


117246 


Combat Sys. Equp. Rm 1 


2-53-1 -C 


4000 


Combat Sys. Equp. Rm 2 


2-126-2-C 


23285 


Combat Sys. Equp. Rm 3 


1-300-0-C 


4000 


Combat Inform. Center 


1-126-0-C 


12120 


Tomahawk Equp. Room 


2-153-2-C 


1421 


5”54 Loader Drum Room 


1-46-0-C 


41658 


VLS No. 1 


01-78-0-M 


138453 


VLS No. 2 


01-338-0-M 


258853 


Starboard Torpedo Tube 


01-342-1-M 


2430 


Port Torpedo Tube 


01 -342-2 -M 


2430 


Gun Mount 21 


03-115-0-M 


13600 


Gun Mount 22 


02-310-0-M 


13600 
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APPENDIX F. COMBAT SYSTEMS AREA RESPONSES-5698 ELEMENT 

MODEL 



The following plots represent the responses for CSER 1 , 2 , 3 and the SPY radar 
room for the 5698 element model. 
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Figure 69. CSER 1 Velocities (Node 3567) 
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Figure 70. CSER 1 Accelerations (Node 3567) 
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Figure 71. CSER 2 Velocities (Node 3669) 
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Figure 72. CSER 2 Accelerations (Node 3669) 
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Figure 75. SPY Radar Room Velocities (Node 6955) 
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Figure 76. SPY Radar Room Accelerations (Node 6955) 
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